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ABSTRACT 

This  study  analyzes  the  influence  of  internal  heat  generation  on  the  temperature  distribution,  displacement  and 
thennal  stresses  in  a  thick  circular  plate  described  as  0  <  r  <  a,  —h  <  z  <  h  within  a  context  of  uncoupled 
thennoelasticity.  A  direct  transient  thermoelastic  problem  of  circular  plate  because  of  internal  heat  generation  with 
radiation  kind  of  boundary  conditions  is  discussed.  The  third  kind  boundary  condition  is  kept  on  upper  and  lower 
plane  surface  of  the  plate.  Over  outer  curved  surface,  (r  =  a),  temperature  is  kept  at  f(z,  t).  The  controlling  heat 
conduction  equation  was  completed  through  Jinite  Hankel  and  Marchi  Fasulo  integral  transform.  The  outcomes 
were  acquired  in  way  of  BesseTs  function  in  the  form  of  infinite  series.  Mathematical  technique  is  employed  for 
aluminium  (pure)  material  plate  and  the  numerical  solutions  were  discussed  and  represented  graphically. 
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1.  INTRODUCTION 

Thermoelasticity  deals  with  the  effect  of  the  thermal  condition  of  a  solid  due  to  application  of  stresses  and  using  the 
reverse  effect,  which  of  deformation.  Thermoelasticity  has  been  aroused  by  the  variety  of  technology  sciences. 
Modern  advancements  of  nuclear  engineering,  aeronautics  and  astronautics,  as  well  as  a  remarkable  advancement  in 
the  subject  of  aircraft  and  machine  formation  have  lead  to  many  challenges  where  thermal  stresses  acts  as  a  part  of 
major  signihcance.  Unexpected  and  ongoing  improvement  of  attentiveness  within  this  subject  grow  throughout  the 
next  half  of  the  twentieth  century,  where  non-isothermal  issues  of  the  concept  of  elasticity  turns  into  highly 
signibcant. 

To  begin  with,  in  this  area  of  aeronautics,  the  large  velocity  of  contemporary  aircrafts  contribute  to  a 
aerodynamic  heating,  where  create  extreme  thermal  stresses,  decreasing  the  potency  of  aircrafts  formation.  Next,  in 
the  nuclear  area,  the  very  substantial  temperature  also  its  gradients  emerging  within  nuclear  reactor  affect  its  design 
and  operations. 

From  the  engineering  of  jet  and  rocket  motors,  the  maximum  temperatures  related  to  combustion  methods 
would  be  the  source  of  undesirable  thermal  stresses.  Therefore,  the  several  issues  where  the  thermal  stresses  acts  as 
a  signihcant,  at  times  the  critical  function,  happen  in  steam  and  gas  turbines,  in  aviation  formations,  chemical 
technology,  in  nuclear  technology,  in  the  science  of  space  vehicles  and  missiles,  as  well  as  in  ship  building,  in 
which,  curiously  enough  ship  fractures  were  usually  credited  to  thermal  stresses  of  average  strengths.  Number  of 
studies  were  completed  to  forecast  deformation  as  well  as  heat  flow  in  a  continuum  within  the  area  of 
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thermoelasticity  throughout  the  past  couple  of  years.  It  still  requires  high  focus  because  of  its  extensive  applications  in 
diverse  disciplines.  Thermal  stresses  also  act  as  a  very  important  part  within  the  area  of  chemical  technology,  at  the 
technologies  of  space  vehicles  and  missiles.  The  thermoelastic  behaviour  of  solids  symbolizing  the  bases  of  containers  for 
hot  gases  or  liquids,  at  the  bases  for  furnaces,  in  applications  including  turbine  engines,  Hywheels,  gears  etc.,  is  highly 
signihcant. 

Nowacki  [8]  studied  thermal  stresses  in  circular  plate  exposed  to  a  axisymmetric  temperature  distribution  over  the 
top  face  without  temperature  on  the  bottom  face  and  the  circular  edge.  Marchi,  E.  and  Fasulo,  A.  [5]  researched  heat 
conduction  in  business  of  hollow  cylinder  with  radiation.  Roy  Choudhary  [11]  talked  about  the  Quasi-static  thermal 
deflection  of  a  clamped  circular  plate  because  of  ramp-type  heating.  Qian  and  Batra  [10]  analyzed  transient  thermoelastic 
problem  of  thick  functionally  graded  plate.  Gogulwar  and  Deshmukh  [3]  established  thermal  stresses  in  thin  circular  plate 
using  heat  sources.  Additionally  Nasser  M.EI-Maghray  [6]  solved  two-dimensional  issue  of  thick  plate  with  heat  resources 
in  generalized  thermoelasticity.  Deshmukh  K.  C.  et.al.  [  1]  examined  quasi  static  thermal  deHection  of  a  thin  clamped 
circular  plate  because  of  heat  generation.  Gaikwad  and  Ghadle  [2]  analyzed  non-homogeneous  heat  conduction  issue  and 
its  thermal  deHection  as  a  result  of  internal  heat  sources  within  a  thin  hollow  circular  disk.  Khan,  Khalsa  and  Varghese  [4] 
explored  reverse  quasi  static  unsteady  state  thermal  stresses  within  a  thick  circular  plate. 

The  research  deals  with  transient  thermoelastic  issue  of  a  thick  circular  plate  covering  the  gap 
D:  0<  T  <  (2,  —  h  <  Z  ^  h,  because  of  heat  generation  with  radiation  boundary  conditions.  This  finite  Hankel 
transform  and  finite  Marchi  Fasulo  transform  method  was  utilized  to  locate  the  solution  of  the  issue.  The  outcomes  existed 
here  to  be  more  helpful  in  technology  issues  especially,  at  the  conclusion  of  the  condition  of  stresses  in  circular  plate 
constituting  bases  of  containers  for  hot  gases  or  liquids,  in  the  bases  for  furnaces,  etc. 

2.  THE  TRANSFORMATION  AND  ITS  ESSENTIAL  PROPERTY 


2.1  The  finite  Hankel  transform  over  the  variable  rand  its  inversion  described  as  in  Ozi§ik[9]  is 


T  (fm,  z,  t)  =  /ra=0  r '  K0  (fm,  r  ')T  (r ',  z,  t)dr ' 

T  (r,  z,  t)  =  Em=i  Ko  (/m,  r)f  (fm,  z,  t) 

where, 

TZ  { r  \  V2  Jo  (fmr ) 

and  the  eigen  value  fm  are  the  positive  roots  of  the  transcendental  equation 
/o  (f  a)  =  0 

where /n(x)is  Bessel  function  of  the  first  kind  of  order  n. 

The  finite  Hankel  transform  t/dciincd  in  equation  (2.1.1)  satisfies  the  relation 


(2.1.1) 

(2.1.2) 

(2.1.3) 

(2.1.4) 


(2.1.5) 
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2.2  The  finite  Marchi  Fasulo  integral  transform  of  f(z),  —h  <  z  <  h  was  described  in  [3]  as 

F(n)  =  J*hf(z)Pn  ( z)dz  (2.2.1) 

Then 


f(z)  =  Sn=i 


(2.2.2) 


where, 

Pn(z)  =  Qn  cos (anz)  —  Wn  sin (anz) 

Qn  —  an(a i  +  af)  cos (anh)  +  (/?x  —  /?2)sin (anh) 

Wn  =  (/?!  +  /?2)  cos  (anh)  +  (a2  —  a^anSin^anh) 

=  LhhPn2(z)dz  =  h[Qn2  +  M42]+^5^[Q„2  -  Wn2] 


(2.2.3) 


The  eigen  values  an  satisy  the  equation 
[o^acos^ah)  +  sin(ah)]  x  [/?2  cos (ah)  +  a2asin(ah)] 

=  [a2acos(ah)  —  /?2  sin(ah)]  x  [f)t  cos (ah)  —  a^asin^ah)]  (2.2.4) 

a^a^.P^  and  (i2  are  constants. 

The  sum  in  (2.2.2)  must  be  taken  on  n  corresponding  to  the  positive  roots  of  the  equation  (2.2.4) 


The  transform  (2.2.1)  has  the  following  property: 


rh  92/(z) 
J-h  dz2 


Pn(z)dz 


Pn(h) 

ax 


Plf(z)  +  «1 


df(z) 


dz 


z=h 


Pn(~h) 

«2 


f2f(z)  +  a2 


df(z) 


dz 


z=—h 


an2F(n) 


3.  STATEMENT  OE  THE  PROBLEM 


Consider  thick  circular  plate  of  radius  a  and  thickness  2 h  covering  those  gaps 

D:  0  <  r  <  a,  —  h  <  z  <  h,the  material  is  homogenous  and  isotropic.  Initially,  the  circular  plate  is  at  zero 
temperature.  On  upper  plane  surface  (z  =  h),  the  third  kind  boundary  condition  was  kept  at  (r,  t),  whereas  on  lower 
plane  surface  (z  =  —h),  the  third  kind  boundary  condition  was  kept  at  g2  (r,  t),  where  (r,  t)and  g2  (r,  t)  are  known 
functions  of  r  and  t.  On  outer  curved  surface  (r  =  a),  temperature  is  maintained  at  f(z,t).  For  time  t  >  0,  heat  is 
generated  in  the  plate  at  the  rate  Q(r,z,t)  .  The  differential  equation  governing  the  displacement  potential 
function  (I>(r,  z,  t) as  in  Nowacki[8]  is 

V2<D  =  KT  (3.1) 

where,  the  Laplacian  operator  V2is 

,  d2  1  d  d2 

V2=  — =  + - +  — = 

dr 2  r  dr  dz 2 
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and 


K  = 


a 


where  K,  vand  a  are  restraint  coefficient,  Poisson’s  ratio  along  with  linear  coefficient  of  thermal  expansion  of  the 

material  of  the  plate  respectively.  <f>  is  Goodier’s  thermoelastic  displacement  potential  and  T  is  the  temperature  distribution 

of  the  plate  fulfilling  the  differential  equation  as  Noda  [7]  is 

+  QG,z,t)  _  1  dT 

X  k  dt 

(3.2) 

Subject  to  initial  condition 

T(r,z,t) 1 1=0  =  0,  0  <  r  <  a,—h  <  z  <  h 

(3-3) 

The  boundary  conditions  are 

T\r=a  =  f(z,t) ,  —h<z<h,t>  0 

(3.4) 

hfT  +  — 1  =  gt(r,  t),  0  <  r  <  a,  t  >  0 

3z  1  z=h 

(3.5) 

h2T +  ^-\  =  g2(r,t),  0  <  r  <  a,  t>0 

dz\z=~fi 

(3.6) 

where  k  is  thermal  diffusivity  of  material  of  the  plate,  A  is  thermal  conductivity  of  material  of  the  plate, 
Q  is  the  internal  heat  also  hlt  h2  are  heat  transfer  coefficients. 

The  displacements  in  terms  of  O  and  M  as  in  Noda  [7]  are 


3<I>  d2M 

Ur  ~  dr  drdz 

(3.7) 

3<I>  .  \n?  *  /i  d2M 

uz  =  -  +  2(l-v)^M-  — 

(3.8) 

M  must  satisfy  the  equation 

V2V2M  =  0 

(3.9) 

The  components  of  stress  as  in  Noda  [7]  are 


aee 


2  G\ 

[©-«•) 

l+£( 

(3.10) 

2  G 

1  M 

CD 

-i  I  e 

1 

)+i 

KM-}f)} 

(3.11) 

2G\ 

[©-«•) 

l+£( 

;(2-v)V“M-g)} 

(3.12) 

2G\ 

[£+£(< 

(3.13) 

where  G  is  shear  modulus. 

For  traction  free  surface,  stress  function 
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ct0/  —  <7rS  =  0  at  z  —  +h  for  thick  plate. 

Equations  (3.1)  to  (3.13)  constitute  the  mathematical  formulation  of  the  boundary  value  problem 
of  thermoelasticity  under  thought. 

z  =  h 

2h 

z  =  -h 

Figure  1:  The  Geometry  of  the  Problem. 

4.  SOLUTION  OF  THE  PROBLEM 


Applying  rule  described  in  above  [2.1]  to  the  equation  (3.2),  and  using  (2.1.5),  we  have, 


<2(fm,z.Q  _  ldT 


(4.1) 


where  T  denotes  the  fmite  Hankel  transform  of  T  and  fmis  the  Hankel  transform  parameter. 


Additionally,  applying  rule  described  in  [2.2]  to  above  equation  (4.1), 


dT* 

dt 


+ kq2T*  =  Y 


(4.2) 


where. 


q2  =  +  a n 2 


and 


Y  =  k{pnQi)§i  ~  Pn(-h)g2  +  Q  (^,n,t)| 

where,  7*denotes  the  Marchi-Fasulo  transform  of  T  and  n  is  the  Marchi-Fasulo  transform  parameter. 
Solving  equation  (4.2),  which  is  a  lst  order  differential  equation,  we  have, 

t)  =  e-^^Ye^dt'  +  c)  (4.3) 

From  (3.3)  and  (4.3),  C  =  0. 

Hence, 

f*(fm,n,t)=  e-k«2t[j*Yek«2t'dt'}  (4.4) 

Applying  rule  in  (2.2.2)  to  the  equation  (4.4),  we  get. 
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T(fn,z,t)  =  S=i  ^e-^^Ye^dt' 


Again,  applying  rule  in  (2.1.2)  to  the  equation  (4.5),  we  get 

7Yr  7  —  lHy  y  ,0  (tm1')  Pn  (z) 

U  ,  ,  ;  a  LmLn  Jo’(U*)  U 


(4.5) 


(4.6) 


where, 


vp  =  e_k<72£|J  rekq2tdt 'j 

and  the  eigen  values  fm  are  the  positive  roots  of  the  transcendental  equation /0  (fa)  =  0. 

The  expression  in  (4.6)  gives  the  temperature  distribution. 

Let  us  consider  MichelTs  tiinction  M,  fulfilling  condition  (3.9)  as 

M(r,z)  =  (/m;)Pn/2)  (4.7) 

Using  (3.1)  and  (4.6),  we  get  displacement  potential  <D  as 

®  =  TR^2nJor^PnxMS(t)  (4.8) 

a  J 0  VS maJ  An 

where 

R  =  k(j—^ja,  and5(t)=  / 'Pdt 

5.  DETERMINATION  OF  DISPLACEMENT  EUNCTION 

Substituting  equations  (4.7)  and  (4.8)  in  equation  (3.7)  and  (3.8),  we  get 


_  V2  y,  y,  tmlo  fimr  ) 

llr  —  2_tm  2-tn 


V2 

u7  =  — 


Jo  (fma) 

Jo  (<fmr)  (z) 
Jo  (fma)  ^n 


fi^5(t)-!L<5 


5(t)- zmz. 


7o(Ur)  pn"  (z)  ) 
‘  Jo  (M 


+  2(1  —  v) 


Sml ]n(fn 


2  Jo  (f mr )  _|_  1  %mJo  P n(z) 


+  SmSn 


r  7o’(U«)  / 

7o(fmr)  (z) 


l/o(Ua)  U 

Substituting  equations  (4.7)  and  (4.8)  in  equations  (3.10)  to  (3.13),  we  obtain 


(5.1) 


(5.2) 


2V2G 


-1? 


1  y,  y,  tmJo  (fmr)  U(z)  _i_  V  V  /o(fmr)  U(z) 

2j?n  2,n  ,  V  e  i  '  2-i  m  n 


X 


+v 


Jo  (fm“)  ^n 

y  y,  (in  Jo  (tmr)  ,  1  y,  y, 
Zj,nZjn  Jo  (fm°) 


Jo  (fmQ)  U 

(mJo  (fmr)  \  rn  (z) 
Jo’(fm«) 


5(0 


,  y1  y1  7o  (fmr)  ?n  (z) 
^mZjn7o'(fma)  An 

fm  /o  (fmr)  P n  (z) 


-SmSn 


/o  (fm“)  ^n 


(5.3) 
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2\/2G 

Oqq  —  ——  X  \ 


-R 


y  y*  *fm  7o  (^mr)  ^n(z)  ■  y1  y1  7o(fmr)  ^n(z) 

2-tm2-tn  ,  Vf  '  2-tm2-tn 


+v 


7o(fma)  +1  m  nl0iUa)  A„ 

fm  7o  (fmr)  ,  ly  y  fm7o  (fmr)  \^n(z) 

'  Zjm  Zjn 


5(0 


7o  (fm<0 

+ 


7o(fmr)  P n(z) 


-iy  E 

Z-im  Zjn 


7o  (fma)  Zn 
^mJo  (^mr)  ^n  (z) 

7o  (fma)  Zn 


(5.4) 


2\/2G 


-/? 


y  y  *fm  7o  (^mr)  ^n(z)  ■  1  y  y  fm7o  (fmr)  ^n(z) 

2-tm2-tn  ,  v?  Zjm  Z<n 


S(t) 


x 


+  (2  -  v) 


7o  (fma)  Zn  r  7o  (fma)  Z 

fm  7o  (fmr)  _i_  £y  y  fm7o  (fmr)  \  ?n  (z) 

'  Zjm  Zjn 


XmXn- 


7o  (fma) 


7o  (fma) 


,  y  y  7o  (fmr)  ^n(z) 

'  Zjm  Zjn  ,  v 


■2m2n 


7o  (fma)  Zn 
7o  (fmr)  (z) 

7o  (fma)  Zn 


(5.5) 


2\/2  G 

OW  = - X 


+  (1  -  v) 


ny  y  %mJo  (fmr)  (z)  y  y  fm  7o  (fmr)^n(z) 

^  2-tm2-tn  ,  uc  _*»  3  *^\L/  2-tm2-tn' 

Jo  ssmaJ  An 


7of(fma)  Zn 


3  2  f  \ 

y  y  £m  7o  (£mr)  ■  1  y  y  £m  7o  (^mr) _ y  y  fm  7o  (fmr)  \  P n(z) 

2-tm2-tn  70'(fma)  r22-tm2-tn 


+  2mS„ 


7o'(fma) 

£m  7o  (£mT)  ^n  (z) 

Jo'(%ma)  Zn 


7oA(fma)  /  Zn 


(5.6) 


6.  NUMERICAL  RESULTS  AND  DISCUSSIONS 


Set/(z;  t)  =  i4(  1  —  e  c)z; 

5i(r,t)  =  52(r,t)  =  B(  1  -  e_t)r, 
Q{r,z,  t)  =  5(r  —  r0)5(z  —  z0)8(t  —  t0) 
where  5  is  the  Dirac-delta  function. 


The  numerical  problem  has  been  completed  for  aluminum  (pure)  plate  with  parameters  a  =  2m,  h  =  lm,  t  =  1  hr. 

Thermal  diffusivity  /c  =  84.18  X  10_6?n2s_1 

Thermal  conductivity /l  =  204.2W /mk 

Density  p  =  2707  kg/m3 

Specibc  heat  Cp  =  896  J/kgK 

Poisson  ratio  v  =  0.35 

Coefficient  of  linear  thermal  expansion  a  =  22.2  X  10-6  //V 
Lame  constant  g  =  26.67 

Also,  =  1.2024,  =  2.7601,^  =  4.3269,f4  =  5.8958,fs  =  7.4655,f6  =  9.0356  are  the  positive 

roots  of 


the  transcendental  equation /0  (2f)  =  0  as  in  [9] 
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r(m) 


Figure  2:  Temperature  Distribution. 


Figure  4:  Radial  Stress  Function. 
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Figure  5:  Angular  Stress  Function. 


Figure  6:  Axial  Stress  Function. 


Figure  7:  Stress  Function. 


From  figure  2,  it  is  found  that,  temperature  turns  maximum  at  the  centre  and  reduces  to  the  external  circular  edge. 
As  time  passes,  the  temperature  of  plate  rises  because  of  internal  heat  source. 

From  figure3,  the  displacement  increases  to  the  external  circular  edge. 

From  figure  4,  the  radial  stress  turns  nothing  at  the  external  circular  boundary  of  plate. 

The  radial  stress  is  compressive. 

From  figure  5,  the  angular  stress  go  on  increasing  from  the  central  part  to  the  external  circular  edge. 
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It  is  zero  at  r  =  0.6.  The  angular  stress  evolves  as  compressive  stresses  near  the  center  and  tensile  stresses 
near  the  external  circular  edge  of  the  plate. 

From  figure  6,  the  axial  stress  varies  nonuniformly  with  radial  direction  and  evolves  compressive  stresses 
in  radial  direction. 

From  figure  7,  stress  varies  with  radial  direction  because  of  effect  of  internal  heat  source.  It  takes  value  zero 
at  the  center. 

7.  CONCLUSIONS 

From  the  research  paper,  the  temperature  distribution,  displacement  and  thermal  stresses  of  a  thick  circular  plate  were 
analyzed  with  identihed  boundary  conditions.  The  signihcant  role  played  by  internal  heat  generation  has  been  represented 
by  Finite  integral  transform  methods.  The  outcomes  were  acquired  in  terms  of  BesseTs  timction  in  the  form  of  series. 
Because  of  internal  heat,  the  radial  stress  develops  as  compressive  stresses  and  the  angular  stresses  evolves  as  compressive 
stresses  near  the  center  and  tensile  stresses  near  the  external  circular  edge  of  the  plate.  The  rate  of  detlection  is  proportional 
to  time  and  increases  as  time  passes.  Substantial  enhancement  in  the  temperature  occurs  due  to  internal  heat. 

Any  specific  scenarios  of  particular  interest  can  be  defined  through  allocating  correct  values  to  the  parameters  and 
functions  in  expressions.  The  outcomes  which  were  acquired  can  be  helpful  to  the  outline  of  structure  or  machines  in 
engineering  applications. 
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